Phase transition in a static granular system 
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We find that a column of glass beads exhibits a well-defined transition between two phases that 
differ in their resistance to shear. Pulses of fluidization are used to prepare static states with well- 
defined particle volume fractions in the range 0.57-0.63. The resistance to shear is determined by 
slowly inserting a rod into the column of beads. The transition occurs at (j) — 0.60 for a range of 
speeds of the rod. 

PACS numbers: 45.70.Cc, 81.05.Rm, 45.70.-n 



A static assembly of granules, for instance sand in a 
rigid container, responds differentlyto shear when packed 
loosely from when packed tightly Q . It is natural to en- 
quire whether these two states are smoothly connected as 
volume fraction varies, or, as with assemblies of particles 
in thermal equilibrium, are such states sharply separated 
by one or more phase transitions. We use recent advances 
in controlling the preparation of granular assemblies to 
show that the latter holds. 

An old magic trick is based on the qualitative dif- 
ference in the resistance to shear of loosely packed and 
tightly packed particles: When a pot with a narrow neck 
is loosely filled with grains, a rod is easily inserted and 
withdrawn. The rod is then inserted and the grains are 
shaken or otherwise agitated to a denser state, where- 
upon the whole apparatus can be lifted by the rod and 
spun about the performer's head 

The existence of distinct phases in granular matter has 
been widely discussed, but a sharp distinction between 
the two phases has remained elusive 0, 0, Eli the dis- 
tinction being hampered by the difficulty in preparing a 
well-defined initial state [6j. The effort to overcome this 
was advanced significantly by Nowak et al. , who used 
a mechanical tapping protocol to obtain well-defined vol- 
ume fractions in the range 0.628 - 0.658. Recently 
Schroter et al. @ showed that a protocol based on ex- 
panding the granular medium by pulses of fiuid from be- 
low could be used to prepare a column of grains with (p 
defined to within 0.1%. Using this technique, we prepare 
granular samples in the range 0.571 < (j) < 0.633. 

Experiment. — To measure the force of resistance to the 
insertion of a rod into a granular sample we use an appa- 
ratus similar to that in Those previous studies 
focused on the influence of geometrical factors such as the 
size of particles, rod, and vessel, and on how the penetra- 
tion force increased when the rod approached the bottom 
boundary. Those experiments were performed at a single 
volume fraction, (j) = 0.59 |^llo|. 

Our measurements are performed with a home built 
granular penetrometer: a translation stage (driven by a 
stepper motor with a step size 2.5 /im) moves a stainless 



steel rod (diameter 6.3 mm and fiat head) downwards 
into a granular sample. The force needed for penetration 
is measured with a load cell with a full range of 10 N 
(Honeywell, Model 31). The sample consists of soda lime 
glass beads from Cataphote with a diameter of 265 ± 
15 /im and a density of 2.484 ± 0.002 g/cm'^ (measured 
with a Micromeritics gas pycnometer AccuPhys 1330). 
The beads are contained in a water-fluidized bed where 
flow pulses of different flow rates allow us to select a 
volume fraction for the static sedimented bed !^ . (If air 
rather than water is used to fluidize a bed, it is difficult to 
obtain low enough volume fraction to see the transition 
[l^.l The beads are fluidized inside a square bore glass 
tube (39.9 x 39.9 mm^). The ratio of inner tube size 
to rod diameter is 6.3, larger than the value five that 
found to be sufficient so that the influence of the vessel 
walls was negligible. Flow pulses are generated using a 
digital gear pump (Barnant Co., model no. 75211). The 
volume fraction is determined from measurement of the 
bed height; one pixel in the digital images corresponds 
to a change of only 0.02% in 0. 

Results. — The force on the penetrating rod increases 
monotonically with its depth, as Fig. ^ (a) illustrates. 
Contrary to the observations in , we find that the rate 
of growth is not polynomial (see the double logarithmic 
plot in Fig. ^(b)). Further, at maximum depth the rod 
is 40 mm from the bottom of the container, well beyond 
the distance of 20 mm where boundary effects have been 
found to be measurable for conditions comparable to our 
experiment |^. 

The rate of change of force with volume fraction ex- 
hibits a well-deflned transition, which occurs for — 
0.598 for a rod depth of 60 mm, as Fig.Elillustrates. The 
(j) values corresponding to the transition increase slowly 
with rod depth (cf. inset of Fig. This dependence 
of the transition point on the penetration depth cannot 
correspond to a dependence of (j) on depth because then (j) 
would have to decrease with depth, which is unphysical. 

There is no hysteresis in the transition: measurements 
at different volume fractions can be made in any order, 
without affecting the results. However, the value of (p at 
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FIG. 1: (a) Penetration force as a function of depth for dif- 
ferent volume fractions, (b) Double logarithmic plot of the 
data in (a). These curves are measured with a penetration 
speed of 10 mm/s, but the results are the same for a range of 
penetration speeds. The total sample height is 110 mm (at <j) 
= 0.6). 



the transition depends weakly on the frictional character- 
istics of the grains, which can change slowly with usage 
of the beads Q • 

The penetration force is independent of the speed of 
the rod (cf. Fig. 12). The absence of a dependence on 
penetration speed is in agreement with the results of 'OJ 
at (/) = 0.59. It also shows that the flow of water in- 
duced by the penetrating rod does not alter the results, 
in agreement with [ilj . 

Typically, only about five percent of the force we mea- 
sure is exerted on the sides of the rod, as we see from 
Figure 121 which measures the force on withdrawal. The 
measurements during withdrawal also suggest the phase 
transition, though it is not well defined since the forces 
are so much smaller. 

The transition between distinct phases indicated by 
Fig. 121 should be accompanied by changes in other prop- 
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FIG. 2: The dependence of the penetration force on volume 
fraction changes at = 0.598 at a depth of 60 mm. Data 
represented by • were measured at a penetration speed of 10 
mm/s, V at 5 mm/min and A at 20 mm/min. The transi- 
tion points are determined by the intersection of least-square 
fits for the volume fractions below 0.595 and above 0.6. The 
inset displays the dependence of the transition point on the 
penetration depth. 



erties. Indeed, measurements of the average height of the 
bed as a function of (/) reveal, as Fig.^lshows, a transition 
at the same (j) as in the force measurements. 

Discussion. — A system similar to ours whose phase 
transitions have been studied is a collection of colloidal 
particles prepared at different </> by centrifugation (13j|. 
The colloidal system has been characterized by a freezing 
density 0/ = 0.407 and a melting density 0m = 0.442, 
such that for (j) < (pf the system is a fluid, for (f> > 0™ 
the system is a crystalline solid, and for (j)f < (p < 0™ 
the system is a mixture of the two states. This phase 
behavior (but with 0/ = 0.494 and (pm — 0.545) was 
found earlier in Monte Carlo and molecular dynamics 
simulations for frictionless hard spheres [l^ . 

There are differences between a colloidal system and 
our static granular system. For example, the colloidal 
system has pure solid and fluid phases separated by a 
coexistence region with a sharp transition at each end, 
while the granular system exhibits only one transition in 
our experiment. Another difference is that the transi- 
tion in the colloidal system corresponds to a change of 
symmetry from disordered to crystalline, while the tran- 
sition in the granular system is marked by a change in 
the resistance to shear. However, a recent consideration 
of the transition in the equilibrium hard sphere model 
|l5l | suggests that the symmetry change may be hiding 
the fundamental mechanism in that model, a type of ge- 
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FIG. 3: Forces measured during a full cycle of insertion and 
withdrawal at a volume fraction 4> = 0.602. The speed of the 
rod is 10 mm/min in both directions. The inset shows the 
withrawal force measured at a depth of 70 mm for the same 
experiments as in Fig. |21 
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FIG. 4: Change in the height of the bed relative to the original 
bed height for a rod inserted to a depth of 60 mm. The 
transition occurs at 4> = 0.597. The bed is slightly deformed 
near the rod but not sufficiently to change the average height 
of the bed. 



ometric constraint, which may also under he the phases 
of the granular system. 

The volume fraction at which we observe a transition 
coincides with earher observations of interesting behav- 
ior. In 8] volume fluctuations around a static steady 
state were measured using a series of identical flow pulses 
in a fluidized bed, and those fluctuations were found to 
exhibit a parabohc minimum at volume fractions between 
0.587 and 0.596, depending on the surface roughness of 
the beads. Arguments based on the central limit theo- 
rem showed that this minimum corresponds to a mini- 
mal number of beads being contained in a statistically 
independent region, which is tantamount to a minimum 
in the correlation length. Another study found that if 
a fluidized bed is slowly defluidized, it exhibits a be- 
havior similar to that of a supercooled liquid, with an 
arrest transition at (/> = 0.594 At the arrest tran- 

sition (j) became nearly independent of the flow rate and 
the correlation length went to zero. A similar result has 
been obtained in an analysis of the Delauney tesselation 
of tomograms of sphere packings: at volume fractions 
between 0.58 and 0.60 re-adjustments involving only a 
single sphere become impossible and any dynamics re- 
quires collective and correlated motion of larger sets of 
spheres (TtI. And finally, as indicated by Figure 0] our 
transition takes place at the critical state density of soil 
mechanics Q; at densities below this the granular sam- 
ple collapses under stress, while at densities above this it 
expands under stress. 

In conclusion, we have taken advantage of a method to 
produce beds of granules with well-deflned volume frac- 
tions to search for a phase transition as a function of 



volume fraction. The transition revealed by our measure- 
ments of penetration force and bed height should help in 
the understanding of granular behavior, for instance, in 
avalanches 0| . The existence of a well-defined transition 
between two phases also suggests the appropriateness of a 
statistical approach to understanding the phases of gran- 
ular matter. 
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